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We have studied the invariant mass distributions of the tttt and KK systems for invariant masses 
up to 1.2 GeV from the J/$ <j)mT{KK) decays. The approach exploits the connection between 
these processes and the tttt and KK strange and non-strange scalar form factors by considering the 
<f) meson as a spectator. The calculated scalar form factors are then matched with the ones from 
next-to-leading order chiral perturbation theory, including the calculation of the the KK scalar 
form factors. Final state interactions in the J/^ — » (f)-K-K{KK) processes are taken into account 
as rescattering effects in the system of the two pseudoscalar mesons. A very good agreement with 
the experimental data from DM2 and MARK-III is achieved. Furthermore, making use of SU(3) 
symmetry, the S-wave contribution to the tt^tt" event distribution in the J/"^ ajTr^vr^ reaction 
is also given and the data up to energies of about 0.7 GeV are reproduced. These decays of the J/\l/ 
to a vector and a pair of pseudoscalars turn out to be very sensitive to OZI violating physics which 
we parameterize in terms of a direct OZI violation parameter and the chiral perturbation theory low 
energy constants and Lg. These constants all come out very different from zero, lending further 
credit to the statement that the OZI rule is subjected to large corrections in the scalar 0'*"'" channel. 
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I. INTRODUCTION 



The decays of the J/^ into a meson and Goldstonc boson pair (tttt or KK) can be used to investigate the 
dynamics of the interacting pseudoscalars. In particular, it was argued in Ref. |Q] that these data together with data 
from pion-pion scattering (and from others sources) force the /o(980) to have a pole structure different to the one 
required by a KK molecule B. This interpretation has been challenged, e.g. in the Jiilich meson-exchange model 
where the /o(980) emerges [Q^ as a KK bound state. Furthermore, in this reference the authors are also able to 
reproduce the data associated with the previous J/vf decays within the same formalism than the one employed in Ref. 
0, but making use of their own strong amplitudes. On the other hand, as will be the topic of this investigation, these 
data can be used to study the violation of the Okubo-Zweig-Iizuka (OZI) rule in the scalar (0++) channel. This rule 
is only well founded in the large iVc limit of QCD, with Nc the number of colors, since OZI violating processes are 
described by suppressed non-planar graphs Q . Still, on a purely phenomenological level this rule works astonishingly 
well, with the exception of the scalar channel, as arg ued e.g. in Refs. [|, 0, §. To be more precise, the decay 
J /'^ 4>MM (where MM denotes the pseudoscalar meson pair) is OZI suppressed to leading order, cf. Fig.|^a, but 
has an additional doubly OZI suppressed contribution depicted in Figj^b. In our approach, both these pieces are 
taken into account. In fact, it will turn out that the second contribution can not be neglected if one wants to achieve 
an accurate description of the data. On the other hand, it is mandatory to have a very precise description of the 
final state interaction in the coupled tttt/KK system (as indicated by the shaded blob in Fig.|l|a) before one can ask 
such detailed questions. As can be seen from Fig.|l|, the crucial ingredient in the reaction at hand are the expectation 
values of the scalar-isoscalar condensates in the pion and the kaon, i.e the so-called scalar form factors. These can be 
calculated at low energies in chiral perturbation theory (CHPT), which is the effective field theory of the Standard 
Model. In our case, the dimeson system can have energies up to 2 GeV and we thus employ unitarity constraints to 
get a precise description of these scalar form factors also at higher energies, demanding furthermore matching to the 
CHPT expressions in the low energy domain. Because of this matching procedure, the large Nc suppressed low energy 
constants L\ and Lg of the next -to-leading order effective chiral Lagrangian can be determined in the process we are 
considering. It has been argued before that so far no direct determinations but rather large Nc inspired estimates have 
been done, see e.g. Refs. 0,0, with the exception of more recent work presented in Refs. [7 10,1^]. Nevertheless, 
as we will discuss in much more detail below, a rather definite determination of L4 can be obtained by considering 
0(p6) CHPT resuhs |l|,|l|. 

To be more specific, to address the problem of the final state interactions in the coupled nn-KK system, we make use 
of the results obtained in Ref. [Q. In this paper it was clearly established that the scattering data of the O"*""*" 7 = 0,1 
(/ denotes the total isospin of the dimeson system) sectors up to centre-of-mass energies of 1.2 GeV are a reflection of 
the strong rescattering effects between the lightest pseudoscalars (tttt, KK for 7 = and tti] and KK for 1 = 1). The 
approach was based on Bethe-Salpeter equations using the lowest order CHPT amplitudes |l^Jl^,^ as the driving 
potential. The fact that one can generate the resonance states for those channels via loop physics, i.e. rescattering, is 
a clear signal of the large deviations from OZI rule in the 0++ sector, see also Refs. 0,^- Such a mechanism has been 
advocated since long, for a pedagogic discussion see Ref. [|l^. On the other hand, it is well known that there is an 
on-going controversy concerning the nature of the scalar resonances /o(980) and ao(980). This controversy originates 
from the observation that there are several different models to deal with the / = 0, 1 scalar sector, all of them 
reproducing the scattering data up to some extend, but with different conclusions with respect to the origin of the 
underlying dynamics. In particular, in Refs. [p^-pO[ these resonances are considered of preexisting origin while in Ref. 

they appear as meson-meson resonances originated by a potential. Also in Ref. [pl| it is advocated for the solution 
that the ao(980), /o(980) are exotic resonances, that is, not simply qq, while the preexisting qq scalar nonet should be 
heavier, around 1.4 GeV or so. Other interesting approaches to this problem are the relativistic quark model with an 
instanton induced interaction of the Bonn group ||2^], the Jiilich meson-exchange approach or the use of QCD sum 
rules [Q. With respect to this controversy, the contribution of the work in Ref. ||lj] is very valuable since, at least, 
the infinite series of diagrams there considered should appear in the whole S-wave partial wave amplitudes calculated 
to all orders in CHPT. The conclusions of Ref. where generalized in Ref. |2^. In that paper, the most general 
structure of a partial wave amplitude when the unphysical cuts are neglected was established. In particular, in this 
paper explicit s-channel resonance exchanges were included together with the lowest order CHPT contribution and the 
whole SU{3) connected scalar sector with / = 0, 1/2, 1 was studied. In particular, it was shown that the amplitudes 
of Ref. ||lj] appear as a particular case when removing the explicit tree level resonance contributions. It was observed 
that the lightest 0"*""*" nonet is of dynamical origin, i.e. made up of meson-meson resonances, and is formed by the 
(t(500), k, ao(980) and a strong contribution to the physical /o(980). On the other hand, the preexisting scalar nonet 
would be made up by an octet around 1.4 GeV and a singlet contributing to the physical /o(980) resonance. With 
respect to this last point, as discussed in Ref. p4], the inclusion of a preexisting contribution to the /o(980) was 
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considered in order to be able to reproduce the data on the inelastic tttt —^ KK cross section^ when including also the 
•qri channel. However, if this channel is not considered, one can reproduce the strong scattering data, including also 
the previous experiments on the inelastic tttt KK cross section, without including such preexisting contribution. 
Finally, in Ref. the contribution in the physical region of the unphysical cut contributions were estimated up to 
■y/s 800 MeV to be just a few per cents making use of the results of Ref. which apply below that energy. We 
will use the formalism of Ref. whose partial wave amplitudes have been also tested in many other reactions. As 
pointed out in Ref. to obtain a consistent picture of the scalar sector, one also has to study other reactions in which 
the O"*""*" amplitudes have a possible large influence via final state interactions. In this way one can complement the 
deficient information coming from the direct strong S-wave scattering data and distinguish between available models. 
In Ref. 1^ all the whole set of photon fusion reactions 77 7r°7r°, tt+tt", K^K~ , K^K and tt^t; were reproduced in 
an unified way from threshold up to ^/s « 1.4 GeV making use of Ref. |l4| to take care of the final state interactions. 
The free parameters present in Ref. were fixed by their values from the PDG In Ref. |2^, making use of 

the formalism set up in Ref. to study the still unmeasured — > ^K^K decay, the reactions jn^n'^, 77r+7r~ 
and jn'^rj were predicted. These predictions were nicely confirmed, almost simultaneously, by a recent experiment in 
Novosibirsk pO| ]. In this manuscript, we will consider yet another applications of the strong amplitudes calculated in 
Ref . jl^ b y studying the J/^E" — > Tnr{KK) decays. In this way, the present study together with the whole set of works 
@,|3g|2| offer a unique theoretical approach to the scalar sector able to discuss all these reactions in an unified 
way. This is achieved without including new elements ad hoc for each reaction, because all these processes are related 
by the use of an effective theory description that combines CHPT and unitarity constraints. 

Our manuscript is organized as follows. In section || we develop and justify a simple phenomenological model for 
the transition of the J/\l/ into the (f> meson and a pseudoscalar meson pair. This model is generalized in section IV 



by making use of flavor SU(3) symmetry and then applied to the w case. Section [II is devoted to the calculation 
of the scalar pion and kaon form factors for the non-strange as well as the strange scalar density to one loop. Then 
we employ the unitarization procedure discussed before to obtain a description of these quantities up to energies of 
1.2 GeV. We also perform the matching of these expressions to the one loop CHPT ones to be consistent with the 
constraints of the broken chiral symmetry of QCD at energies below approximately 0.5 GeV. The results are presented 
and discussed in section [V where we also generalize section ||. Our conclusions and outlook are given in section 

II. MODELING J/* <Ptvtt , (f,K'K DECAYS 

We will calculate the S-wave contribution to the invariant mass distributions of the tttt and KK systems in the 
J/^* (f>inr{KK) decays. Taking care of the final state interactions of three particles can be simplified to a large 
extend if only two of the final state particles undergo strong interactions, and the third is merely a spectator. We 
will assume that this is the case in our present problem and we will take the 4> as the spectator. This is certainly 
very well sounded for the ^tttt system since the (j)TT interaction is very weak as required by the OZI rule. On the 
other hand, the situation is not so clear with respect the (f> and the kaons. Nevertheless, at the energies in which the 
kaons become important, above the KK threshold, the experimental mass distribution is completely dominated by 
the /o(980) resonance and this state is a two body effect emerging from the coupled tttt and KK systems, as discussed 
already in the introduction. Since we are only considering a small range of energies above the KK threshold, this 
approximation should be justified. 

We therefore describe the transition from the J/^ to the (f)+2 Goldstone bosons system by an effective Lagrangian 
based on the following phenomenological arguments: 1) The already discussed spectator role of the (f> resonance and 
2) the TTTT and KK invariant event distributions, which will be shown later, seem to be clearly dominated by the 
S-wave contribution, although these experimental data have not yet been subjected to a partial-wave analysis. These 
experimental facts, together with Lorentz invariance, can be easily incorporated in the formalism just by writing the 
interaction vertex of the J/^ resonance with the (j) meson and some scalar source S with vacuum quantum numbers, 
J^'~' — O"^"*" with spin J, parity P and charge conjugation C, as: 

(1) 



^In the last edition of the PDG tables |25j it is argued that, possibly, the previous experiments have a much larger uncertainty 
than previously given in the corresponding publications. 
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with g a real coupling constant. We briefly discuss why other possible structures involving derivatives on the various 
fields should be suppressed. The J/"^ is very heavy and thus can be considered a static source. Derivatives acting 
on the (j) and the scalar source S can be combined to the invariant structure '^^{d^cf)^ — d^(l)^)d'^ S. This leads to a 
vertex of the form e(^e^P0 • (^i + 92), with qi^2 the momenta of the two Goldstone bosons and the momentum of 
the (f) meson. However, due to momentum conservation, we have • (91 + 92) = {M"^ — — s) /2, with s — {qi+ (72)^ 
the total two Goldstone boson energy squared. Due to the large value of the J/^ mass, this combination of momenta 
is essentially constant for the Goldstone boson energies considered here, ^/s < 1.2 GeV. Such terms become more 
important at higher di-pion (kaon) energies. Therefore, we can generically write such type of higher order corrections 
to Eq.(|l|) in the form 

/ ((91 + '?2)^P0• (91 + 92)) , (2) 

where the function /(. . .) essentially only depends on the first argument. Such terms that depend on (91 + 92)^ can 
be derived from an interaction of the type ^^(j)f_id^d'^ S . Such structures lead to a weak s-dependence of the constant 
g and/or the parameter defined below. We have checked that such (weak) energy dependencies do not change 
any of the conclusions obtained when treating g and as energy independent. Another possible higher order term 
of the form ^ ^(p^d^^d" S giving rise to the vertex e(^e^ (91 + 92)^(91 + 92)"^- Such couphngs can also have an S-wave 
contribution, which can be obtained by properly summing over the pertinent polarization vectors. Again, due to the 
large mass of the J/'^ , such terms are only weakly s dependent and can be treated along the lines outlined before. 
More complicated structures can always be brought into some linear combination of the ones just discussed or have no 
S-wave component. These considerations not only show that our ansatz Eq.(|]) is quite sensible in the energy range 
considered here but also that corrections to it can be worked out consistently. 

From the OZI rule, which can also be seen as a result of the large Nc expansion of QCD and the experimental 
absence of any clue indicating a non negligible interaction between the (j) and the pions, one should expect that this 
scalar source S would be simply made of strange quarks, i.e. S ~ ss. However, it is known that the also decays 
into non-strange mesons and furthermore, there are strong arguments to believe that large violations of the OZI rule 
(and of the large Nc limit of QCD) are manifest in the 0++ sector [[7P, p^ , p^ , as discussed in the introduction. As a 
result, we will consider a more general scalar source S, that also has a contribution of the form X^nn, where 

nn = —= {uu + dd) (3) 
v2 

parameterizes the contribution from the non-strange quarks and A^ is just a constant measuring the relative strength 
of this contribution with respect to the strangeness component ss. Already at this point we stress that the choice 
A^ ^ will be justified a posteriori by the results presented below. Therefore, we use 

S — ss + X^fin (4) 

Consequently, it follows from Eqs.(0) and (jij) that the transition matrix element for the process J/^ ^ (f>MM is 
given by: 

T = e(*; p)^e{(j); p'Y (0| {ss + A^ nn) \MM)* (5) 

where |0) is the vacuum state, e(^'; p) is the polarization four-vector of the J/^ resonance with polarization p and 
analogously e(\E'; p') is the polarization four-vector of the (j) resonance, and the * denotes complex conjugation. Note 
that in this equation we are implicitly assuming that the is a spectator as discussed before. The scalar source S 
couples to the two meson system, in which the rescattering (final state interactions) appear. The anatomy of our 
model is depicted in Fig.^ As discussed below, invoking SU(3) symmetry, we will also apply this approach to the 
S-wave contribution of the J/^ — > wtttt decay to further constrain the description of the measured event distributions. 



III. COUPLED CHANNEL PION AND KAON SCALAR FORM FACTORS 

As a consequence of Eq.(H), our problem is reduced to calculate the matrix elements {Q\ss\M M) and (0|nn| MM), 
which correspond to the strange and non-strange isospin zero (/ = 0) scalar form factors. 
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A. Definition of the scalar form factors 



One can define an extended QCD Lagrangian allowing for the presence of external sources. In this way the 
identification of matrix elements of quark currents can be done easily. For instance, a scalar source can be added 
simply as: 



(6) 



where q embodies the three light quarks, u, d and s. The QCD current quark mass term can be obtained from such 
a scalar source by setting, 



E = diag(m„, md, rus) 



(7) 



This is the standard method of treating explicit chiral symmetry breaking in CHPT (or any similar effective field 
theory). Consequently, we can work out the scalar quark-antiquark operators. 



dd 



dCQCD 


_ dCqcD 




druu 




_ dCqcD 


9E22 


dmd 




dCqcD 


0^33 


drus 



(8) 



One can include in the same way as in the QCD Lagrangian the external sources in the effective CHPT Lagrangian |9[| 
simply based on symmetry arguments. In the lowest order chiral effective Lagrangian, C2, the scalar source appears 
in the mass term 



/•mass 
^-2 



(9) 



with / the meson decay constant (in the chiral limit), x = 2i3oE and Bq is a constant not fixed by symmetry. 
This constant parameterizes the strength of the quark-antiquark condensation in the non-perturbative vacuum, 
Bq = \{0\qq\0)\/ f"^ . The trace in flavor space is denoted by (...). The octet of Goldstone bosons is collected in the 
matrix-valued unimodular field U{x), 




(10) 



with 



71-" + 7E^^ 

TT 

K- 



TT+ K+ 



%/6 



18 



(11) 



It is then straightforward to work out the scalar-isoscalar quark-antiquark operators from the effective Lagrangian, 



dd 



ss 



dC2 
'5S22 

'dll33 



-fBo 



-fBo 



-fBo 



1 - 72 ' 



2 6 VS 



K+K 
K°K 



6 V3 



' 2 



(12) 



where the ellipsis denotes terms of higher order in the meson fields not needed here. From the last of these equations 
one concludes that the strangeness component of the pion should be small since it only comes in at higher orders. 
From this representation of the scalar operators, one can deduce the pertinent expressions for the scalar form factors. 
For our purposes, it is sufficient to consider pure isospin zero (/ = 0) states formed from a pion or kaon-anti-kaon 
pair, i.e. 
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= -^\K"^K- + . (13) 

v2 

Note the extra factor 1/V2 in the definition of the / = Itttt) state, it is introduced to take care that in the isospin 
basis states the pions behaves as identical particles. Combining this with Eq.(^2|), one can easily calculate the lowest 
order (tree level) CHPT results (remember the normalization of the non-strange quark operator given in Eq.(|^)), 

(0|nn|TTTT 
lS>\hn\KK) = Bo , 
(0|ss|7T7t) = , 

{0\ss\KK)^V2Bo. (14) 

As anticipated, to leading order the two-pion system has no strangeness component. To all orders, these matrix 
elements are given in terms of four scalar form factors,^ 

(0|nn|TTTT) = V2BqPI{s) , 
{0\nn\KK) = V2BqT^{s) , 
(0|ss|tttt) = V2BoTl{s) , 
{0\ss\KK) = %/2Bor^(s) . (15) 

Here, the following notation is employed. The superscript s/n refers to the strange/non-strange quark operator 
whereas the subscript 1, 2 denotes pions and kaons, respectively. In the following we will remove from Eqs.(0[5l) the 
overall factor v^Bq, since the experimental data on the J/^ — > (pMM decays are not normalized. 



B. Next-to-leading order pion and kaon scalar form factors 

The pion scalar form factors r"(s) and Ff (s) were calculated in Ref. up to one loop in CHPT. Since they were 
not explicitly given in Ref. 1^, we give here the pertinent expressions: 

rUs) = |l + - + ^ {2Ll - LI) + 8{2Ll L^ ^^K + ^rnl ^ ^^^^ ^ 2 ~ 
P V3 



T{{s) = {2LI - Ll)^^^ + i=/(.) , (16) 



with /(s) and f{s) given by 



2s-ml - , , s - , , ml - , , 4s f,^ 1 /., ml 



/(.) = r_pL J„,,,) _ _r_J,,,(,,) _ ^ J„„) + ^ |i; _ _^ (^41„g - log + 3 

/(.) = I + f .7„„(.) + - 5^ (log f + 1) } , m 



and Jpp{s) {P = TT, K, rf) is the standard meson loop function |^ 

(2 + ap(s)log- 
Idtt^ \ (jp{s) + 1 



Jpp{s) ^ -L- ("2 + ap{s) log ^^i^i ^ ] ) , (18) 



■^We remark that more commonly the definition of these form factors includes the pertinent quark masses, such that e.g. the 
non-strange scalar form factor of the pion is defined via {0|rn {uu + rfd)|7r7r) = r.^(s). For our later discussion, the overall 
normalization does not play a role but should be kept in mind. 
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and ^ is the scale of dimensional regularization. The quantities /ip in Eq.(|T6|) are given by 



Hp 



nip nip 



(19) 



The scalar kaon form factors at next-to-leading order in CHPT are not given explicitly in the literature. We fill here 
this gap by performing such a calculation. This implies calculating the diagrams shown in Fig.|3, which comprise the 
lowest order CHPT result, Fig.^, already derived in section III A, the tadpole contribution, Fig.|b, and the unitarity 
corrections, Fig.||c;. The vertices for these diagrams come from the lowest order CHPT Lagrangian. We note that wave 
function renormalization diagrams are not depicted in this figure. Finally, in Fig.^, the local contribution coming 
from the 0{p*) CHPT Lagrangian is depicted. These terms are parameterized in terms of the scale-dependent, 
renormalized low energy constants L\{h) (in our case i = 4, 5, 6, 8). Evaluating these diagrams leads to 



9s - ^ml 



3s 



36/2 
4L 



8K 



K 



P 



16Z/g / . 2 2 \ ^ 



(20) 



As a test of our calculations we have checked that the infinities, associated with the wave function renormalization 
contributions and the loops in Fig.^ and ^ arc properly absorbed by the infinite parts of the pertinent low energy 
constants and thus the expressions given in Eqs.(^2) s-^e finite. 



C. Unitarity requirements 

We now discuss the constraints that unitarity imposes on the scalar form factors. Of course, at low energies, one 
can simply work with CHPT and treat unitarity in a perturbative fashion. Here, however, we are interested also at 
energies of the order of 1 GeV, which requires some resummation technique and also the channel coupling between 
the TTTT and the KK systems has to be taken into account. This has been elaborated in big detail in Ref. |14| and 
we present here the formalism necessary to discuss the scalar form factors. For convenience, we employ the matrix 
notation already introduced in the previous subsection, i.e. pions are labeled by the index 1 and kaons by the index 
2. ^From Ref. [ |T4| , we have the following expression for the T-matrix for meson-meson scattering, 



r(s) = [/ + K{s) ■ g{s)] ■ K{s) 



(21) 



where s denotes the centre-of-mass energy squared and K{s) can be obtained from the lowest order CHPT Lagrangian, 

s - ml/2 



K{s)i 



P 



3s 



(22) 



We remark that because of time reversal, both K{s) and T(s) are symmetric functions, so that K{s)2i — K{s)i2 and 
similarly for T{s). The matrix g{s) is also diagonal and given by |32f| 



1 



167r2 



cri(s)log- 



a.(s)./l + -^ + l 



a,(s)./l + 



-2 log 




z = l,2. 



(23) 
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where ai{s) = ^1 — Amf/s, /tt ~ 93 MeV is the weak pion decay constant, rrii are the masses of the pions (mi = 138 
MeV) and kaons (m2 — 495.7 MeV) and gmax = 0-9 GeV is a cut-off in three-momentum space. On the other hand, 
g{s)i can also be calculated in dimensional regularization, using the standard MS — 1 scheme employed in CHPT 

= tAi f-1 + log 4 + ^^(^) = -•^^(^) ' (24) 

where the last equality follows from the definition of the renormalized two-meson loop function p] 



1 

167r^ y ' " 



Ms)-7^\^ + ^^S^]=Ju{s) ■ (25) 



Note that we have changed the subscript ^^pp" appearing in Eq.([l^) into "ii" to conform with our matrix notation. 
By expanding Eq.(p3|) in terms of mi/gmax one can easily see, as discussed in appendix 2 of Ref. that the 

differences between g{s)i in Eq.(^3|) and Eq.(|2j) are of higher order in the chiral expansion, i.e. of order 0{m^/q^^^), 
for the following value of the scale fi, 

H = — 1=- ^ 1-2 gmax • (26) 

For energies above the threshold of the state i, unitarity implies the following relation between form factors and the 
/ = T-matrix: 

Im r.(.) = ^r,(.)^0(. - 4mJ)(7;S-wave(^))* (27) 



with Pi{s) = — mf the modulus of the cm. three-momentum of the state «, and the strong amplitudes are 

projected on the S-wave. In the former equation we have suppressed the superscript "n" or "s", appearing in Eq. (p^ , 
since the previous equation applies to any of them. Finally, in what follows, we will also remove the superscript "S- 
wave" with the understanding that any partial wave is projected onto the S-wave. Taking now the complex conjugate 
on the right-hand-side of Eq.^^) and using the fact that the T-matrix is symmetric, we can rewrite Eq.(p7|) in matrix 
notation as: 

Imr(s)=T(.).g^.r(s) (28) 

where 

Substituting in the previous equation Im V{s) by (r(s) — r(s)*)/(2«) and T{s) by its expression given in Eq.(|2l]), one 
has: 

T{s) = |/ + [/ + K{s) ■ 9{s)]-' ■ K{s) ■ • T{sr 

= [I + K{s) ■ g{sr' ■ |/ + K{s) ■ g{s) + K{s) ■ • r(.)* . (30) 



Taking into account that the A'(s)-matrix, Eq. (|22|), is real and that 



since 



■''■>+'S^ (31) 



we can write Eq.(|30|) as: 
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[/ + K{s) ■ g{s)] ■ r{s) = [I + K{s) ■ gis)*] ■ r(s)* . (33) 

This tells us that the quantity [/ + K{s) ■ g{s)] ■ T{s) has no cuts since the only one which appears in g(s) and r(s), 
the right or unitarity cut, is removed. Therefore, we can express T{s) as: 

T{s)^[I + K{s)-g{s)r' -Ris) (34) 

with R{s) being a vector of functions free of any singularity. We remark that this procedure of taking into account 
the final state interactions is based on the work presented in Ref. | |3^ . In the following, we will fix R{s) by requiring 
the matching of Eq.(^) to the next-to- leading order (one loop) CHPT tttt and KK scalar form factors. These are 
calculated in the next subsection. 

It is worth to stress that Eq.(^4|), given in terms of a vector of functions R{s) without any cut, can be applied to 
any K-matrix without unphysical cut contributions, as the one derived in ref. p4| . The use of the strong amplitudes 
calculated from this reference is appealing for several reasons: 1) Because of their simplicity, 2) they have been already 
successfully used to describe many two meson production processes, as discussed in the introduction, and 3) higher 
order corrections to the kernel used in ref. |l4[| are not necessary to match with the next-to-leading order CHPT scalar 
form factors. In fact, in ref. p3 one can find a detailed comparison between the approach of ref. and the more 
general ones described in rcfs. [P4| and ||3^. The main conclusion is that, apart from the detail of including (or not) 
the rjr] channel as already discussed, the unitarity corrections coming from the rescattering of the lowest order CHPT 
kernel completely dominate the strong S-wave / = scattering amplitudes up to about 1.2 GeV. Thus, one would 
not expect relevant departures from the use of the strong amplitudes from ref. |jl^ or from refs. or |2^. In fact, 
all these approaches give rise to very similar pole positions for the /o(980) and a mesons. For higher energies new 
effects have to be taken into account as e.g. the contributions from a pre-existing octet of scalar resonances around 
1.4 GeV 00 and the increasingly important role played by multiparticle states, basically the An intermediate state. 
In addition, one has to deal with more relevant interaction vertices between the various fields than those given in 
Eq. (|l]) as discussed in section 

D. Matching with chiral perturbation theory 

The general expression for the scalar form factors given in Eq. (jsj) can be further constrained by matching it to the 
one loop CHPT expression given in Eqs.(^, |2^). This ensures that for energies where CHPT is applicable, these form 
factors fulfill all requirements given by chiral symmetry and the underlying power counting. This matching procedure 
essentially fixes the vector R{s). We remark that since in our unitarization procedure we are not considering the ijij 
channel we thus can not reproduce the chiral logarithms associated with this channel. Therefore, we will only consider 
the contribution form this channel to the value of the form factors at s = and we will not include any s dependence. 
This approximation should not induce any sizeable theoretical error because the influence of the 7777 channel was found 
to be significant only above its threshold [p4|j3^ ] as already discussed in the introduction (when comparing the results 
of Ref. and Ref. 

We only discuss in detail the matching for the form factor r|(s). The procedure for the other form factors is completely 
analogous and we thus only give the final results for them. From Eq.(Q) one has: 

r^(s) = [/ + Kis) ■ gis)r' ■ R^{s) - [/ - Kis) ■ 5(5)] • R%s) + 0{p') (35) 

where the superscript "s" in R'^{s) indicates that we are considering the ss form factor. From the former equation 
and Eq.(p^) one has that R'^{s)i = 0{p'^) and that i?*(s)2 = 1 + 0{p'^). Hence, we can recast Eq. (p5[) as: 



rish = R'ish - K{sh, gish + 0(/) = R'is)2 + ^ JW{s) + 0(/) (36) 



at the regularization scale /i = 1.2g„iax- Comparing this result with the one given in Eq.(|2C|) leads to 
R'{sh = 1 + ^(-s - 4mi.) + ^(. - 4mi. - m^) + + ^(4m^^ + + \ 



a + log^), (37) 



3 



Preexisting means here that these resonances with a mass around 1.4 GeV are as "elementary" as the basic fields vr, K or r). 
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using the Gell-Mann-Okubo relation = 4m|./3 — »7iJ/3, the deviations from it being of higher order for our 
purpose. Proceeding in an analogous way for the other form factors one concludes: 



/3J A{Ll + 2Ll) ^ m{2Ll-L,) 
2] p p 



8(2£g - LI 
P 



■{2m 



K 



3m, 



(6m 



m 



K 



72 7r2 f 



■(1 + log^) 



i?^(.)i=V3^(.-2m^) + i^m^ 



(38) 



Notice that R'^{s)i is subleading in large Nc, i.e. of 0{N^^), while the other quantities in Eqs.(37[3^) are of order OCV) 
in this counting. This is expected since the produ ction of pions from an ss source is subleading in large Nc QCD [||. 
We also see, as already stressed in section [IIA that R'^{s)i is 0{p^) in the chiral counting. Once the functions 
i?"''*(s) have been determined, the final expressions for the form factors are obtained by making use of Eq.(Q |2|, 
p3| ). Finally, one has to take into account that, when using Eqs.(|37|, ^8| ), the regularization scale is fi — 1.2gmax — 1-08 
GeV. Therefore, we have to run the low energy constants L\{^) to this scale from the usual ones /i ~ rriri or fj, — nip, 
with rrirj, rup the mass of the r], p meson, respectively, by using the appropriate /3-functions given in Ref. M. 



IV. RESULTS 

We will first discuss the results for the J/^E* (j)TnT{KK) decays and then we will also consider to some extent the 
J/^ (jjTTTT decay. To be more specific, we consider the S-wave contribution to these decay modes. 



A. The <A-meson case 



Considering the phase space of three particles [ p5| we can write the unpolarized event distribution for the J/^ 
bn^TT^ {K^ K~) reactions as: 



dN{W)^ 
dW 



(2^)312 m2/^ 



|rns) + A^rr(s)|= 



1 ^ o o 



(39) 



and i = 2 to the K^K^ system, in order. Furthermore, W is the total energy in 
the cm. of the two pions or kaons, Gi is basically a Clebsch-Gordan coefhcient equal to 4/3 for pions and 1/2 for 
kaons, respectively, and a normalization constant depending on the experiment, in our case DM2 or MARK- 
III [^^. In comparing with the experimental data, we will average Eq.(|3^) over the width of the bin (as given by 
the corresponding experiment). As discussed in section HID , our calculated form factors depend on the CHPT low 
energy constants L\, L^, Lg 



and ig. From these, only Lg and Lg are relatively well determined. Their most recent 



values, given in Ref. [g2[ from an 0{p ) CHPT analysis of the Kgn form factors, are: 

10^ Ll{Mp) = 0.65 ± 0.12 , l{y^Ll{Mp) = 0.48 ± 0.18. 
At the scale = 1.2(7max — 1-08 GeV they are: 

10^ L^(1.08 GeV) = -0.15 ± 0.12 , 10^ Lg(1.08 GeV) = 0.26 ± 0.18 



(40) 



(41) 
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1^ as given in Eg. (p^) . Thus, although at present no precise determination of this 
ler strong constraints on its value can be obtained by combining the determination 



On the other hand, L4 and Lg are only poorly known and their present values Q can be considered as stemmiiig more 
from an estimation of their order of magnitude than from a truly phenomenological fit. According to Ref. [y, it is 
estimated that for a regularization scale fi between 0.5 and 1.0 GeV one should have 10'^ L4 ~ ±0.5 and 10'^ Lq ~ ±0.3. 
A more recent determination gives 10^ L\{mp) = —0.3 ± 0.5 and 10'^ LQ{mp) = —0.2 ± 0.3 so that at a scale of 
1.08 GeV one has 

103l:;(1.08 GeV) = -0.57 ±0.5 , 10^ ig(1.08 GeV) = -0.36 ± 0.3 . (42) 

Again, this estimate relies on OZI (large Nc) arguments. To be more precise, one sets L\ g to zero at the scale 
/i = m^, which is of course somewhat arbitrary. One can also make use of information about the low energy coupling 
constant £4 coming from two flavor CHPT by means of the relation £\{^j,) — 8L\{ii) + 4Lg(/x) — vk/'^ + 0{p^) with 
vk = [ln(Af|-//x^) + l]/327r^. The low energy couphng constant has been determined at 0{p'^) in ref. |^ with the 
result 10^ £4 = 1.2 ± 6 at the scale ^ = 1.08 GeV. However, in ref. Q making use of the analytically deduced pion 
scalar form factor uu + dd up to and including 0{p^), they update the previous value and give the improved result 
lO'^ £4 = 1.8 ± 1.9 at the same scale. Although the central value from both determinations is very similar the error is 
much smaller in the second case. With the value for given in Eq. (|4l|), we obtain 10'^L4(1.08 GeV) = 0.1 ± 0.7 
when using £\ from ref. Q and 

10^ L^(l. 08 GeV) = 0.19 ± 0.25 (43) 

in the case of 0{p^) SU(2) CHPT jl^. The value obtained in this way for L\ is also consistent with zero, but on 
the positive side, in stark contrast to the value given in Eq.(^) whose central value is very far from the more precise 
determination coming from refs. 
low energy constant is available, rat 
of the low energy constants making use of 0(/) CHPT, both in its SU(2) [ij and SU(3) ||12| forms. 

More recently, some constraints on the couplings Kand Lg have been reported ||ic| , [lll,|7|] . The determination of L\ 
relies on a comparison of the CHPT series at next- |£0| and next-to-next-leading order 111]] with a phenomenologically 
determined scalar form factor. In ref. the value 10'^L4(1.08 GeV) ~ 0.14 is given without errors and the band 
of values -0.12 < 10^ L!;(1.08 GeV) < 0.04 is reported in ref. Q. We note that the second set of values for LI 0] 
is in the lower limit of the value for L\ given in Eq.(^). Nevertheless, the determination of ig ||l^,|ll| is not so well 
sounded due to strong simplifying working assumptions when computing phenomenologically the quark correlator 
{uu + dd)ss @,|lll- It is also stated in ref. |lOf| that the positivity of the fermionic measure gives rise to a lower bound 
for Xg, 10'^Lg(1.08 GeV) > 0.03. However, this bound is somewhat arbitrary since the only necessary requirement 
to make use of the positivity of the fermionic measure is that all the three light quark masses have to be equal. In 
ref. they were set equal to the strange quark mass, but they could as well have been taken on another value. In 
fact, for the average light quark masses {rriu + md)/2, the corresponding lower bound is: 10'^Lg(1.08 GeV) > —0.75. 
The bound based on using the strange quark mass can only be maintained if one assumes the next -to-leading order 
corrections to be of canonical size, cM^/(47r/Tr)^, with c a number of order one.^ 

While the data of DM2 |^l| have been published, this is not the case for the data of MARK-III . On the other 
hand, both experiments, see FigsJ^ and |, are compatible for the tt+tt^ event distribution for 25 and 10 MeV bins^. 
However, this is not the case for the KK event distribution, see Fig.|^. Consequently, we will only consider for our 
fits the data from DM2 |^ for the tt+tt" event distributions in the J/^ 07r+7r~ decay, both for the 10 and the 
25 MeV bins. In fltting the data for the J/^E* ^tt+tt" distribution from DM2 we will fix Lg and Lg at the values 
given in Eq.(|4^). On the other hand, L\ and Lg will be taken as free parameters. In this way, our expression for the 
event distribution of the pions and kaons will have four free parameters: C^, A^, L4 and Lg. However, there is still 
too much freedom in fitting the data with this set of free parameters since the global normalization constant can 
only be determined within a large range and with a sizeable uncertainty. We will further restrict our fit by requiring 
that we can also describe the pion event distribution in the J —> wtt^tt" decay, at least in the low energy region 
where the S-wave contribution, the one which we are considering here, is dominant. This extension of our model to 
the Lo case is discussed in the next subsection. Imposing this requirement, is fixed,QC^ = (16 ± 3) MeV~^, and the 
fitted values for A^, LJ(1.08 GeV) and Lg(1.08 GeV) are: 



*We are grateful to Bachir Moussallam for a claryfing discussion on this topic. 
^The data for the 10 MeV bins of MARK-III has been taken from Ref. 

^There is approximately a factor 1.11 between the global normalization constant required for the DM2 data with respect the 
one required for MARK-III. 
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= 0.17 ±0.06 , 10^ L^(l. 08 GeV) = 0.44 ±0.11 , 10^ L^(1.08 GeV) = -0.38 ± 0.06 , (44) 

with a x^/dof = 0.92. Clearly ^ 0, in contradiction with the OZI rule. Furthermore, the pion event distribution 
turns out to be very sensitive to the large subleading low energy constants L\ and Lg. The theoretical uncertainties 
given in Eq.(Q) are obtained in the following way. We have allowed for a relative change of 20% in the global 
normalization constant when considering the data with the and also the ones with the omega in the final state. 
We consider this estimate of the error in as conservative, since the ensuing deviation from the wtt+tt" data for 
such changes in is larger than the uncertainty in the omega data by assuming a Poisson distribution. On the 
other hand, we have also allowed an uncertainty of ±0.1 GeV in the determination of gmax from ref. |l^ and then we 
calculated the band of values for A^, L\ and Lg. All these sources of uncertainty are added in quadrature together 
with the statistical error given by the fitting procedure when using the central values for and gmax- 

It is instructive to compare the values for the low energy constants found here, e.g. Eqs.(Q), with the ones given 
in Eq.(^2|). While ig agrees perfectly within error bars, the sign of L\ is changed. Stated differently, if we evaluate 
from Eq. (^^ 10'^ L\ at the rho mass, we find a value of 0.71, which is sizeably larger in magnitude than the central 
value given in Ref. |Q. Also, it is larger than the positive value deduced from SU(2) information given in Eq.(^3|), 
at the scale /i — nip one has for this case 10^ L\ — 0.46 ± 0.25, although both values are consistent within errors. We 
reiterate that using large arguments, one would expect L\ to be zero (at a scale somewhere in the resonance region). 
Therefore, our increased value and also the one from Eq. (^3|) , clearly signals OZI violation. Quite differently, our value 
for lO^Lg(mp) = —0.22 is completely consistent with the previous determination [Q. However, it does not fulfill 
the positivity constraint lO^Lg(mp) > 0.20 but it fulfills the other reasonable 'positivity' bound previously 

discussed lO'^Lg(mp) > —0.59. In fact, our value for 10'^ ig lies in an natural intermediate region between both 
extreme lower bound. Nevertheless, we have also performed a series of fits enforcing the former constraint. We can fit 
the (/) data, but on the expense of very large values for A^ and L\. Furthermore, it is not possible to simultaneously 
get a description of the w decay data. We think that further study is needed in order to apply the positivity of the 
Dirac measure and also, we should stress that the LEG Lg is plagued by the Kaplan-Manohar ambiguity It is 
important to point out that one can criticize our determinations of the low energy constants for two reasons. First, 
our model for the J/^* decay with the meson as a spectator is fairly simple, one could e.g. write down higher order 
transition operators which would complicate the analysis. Given, however, the fact that we can precisely reproduce 
the data both for the and the w resonances, it is not obvious a priori that such a modified ansatz would lead to very 
different results. Second, the use of unitarity to determine the scalar form factors beyond one loop accuracy induces 
some inevitable model dependence. To overcome this, one could think of doing a pure CHPT analysis on the left wing 
of the scalar resonance. We believe, however, that the present data in this energy region are not precise enough for 
an accurate determination of the LECs. Independently of these reservations, our analysis clearly underlines that the 
OZI rule is strongly violated in the scalar 0++ sector, as indicated e.g. by the large positive value of the L EG L \{mo) 
and also by the non- vanishing value of A^. With respect the latter point, see also the footnote in section 



IV B 



It is also worth to indicate that in ref. ||3^, making use of the Inverse Amplitude Method (lAM) |Q with complete 
next-to-leading order CHPT strong amplitudes, a fit to the / = and 2 S-wave and 7=1 P-wave tttt and KK partial 
wave amplitudes was done in terms of the low energy coupling constants L\, L2, L^, L\, Lf and 2Lg ± Lg. This study 
has in common with the present one that a complete matching to the relevant next-to-leading order CHPT results was 
given and at the same time fully unitarity amplitudes were derived. The experimental data was very well reproduced 
up to energies around 1.2 GeV giving rise to the presence of the resonances p and /q(980). The values obtained for 
Lg and 2Lg + Lg, within errors, are consistent with those recently obtained in ref. |lj]. This implies agreement of 
the results of that reference with our choice for the values of L^and Lg |l^] and our presently determined value for 
Lg. The main difference between the set of values given in ref. [Q and those in ref. ||l2[ corresponds to the value of 
L2- While in the former case L2 — 2Li as required by Vector Meson Dominance (VMD) F] this relation is only 
fulfilled within errors by the values given in ref. |l|. Nevertheless, in ref. || 10^ L^(TOp) = 0.2 ±0.1. This value, 
although on the positive side, is incompatible with our present one, lO'^Kfmp) — 0.71 ± 0.11, and, within errors, is 



compatible with the rest of analyses presented in this section except for ]ll| ]. Summarizing, for Lg there is a rather 
good agreement between our present study and refs. |0, jl^, |Q in disagreement with the finding of ref. ||l^,|ll|. On 
the other hand, for L4 our present analysis is compatible only with that value of L\ determined from 0{^) CHPT 
which also find quite a sizeable central value at /it = Mp, around 0.5 x 10"'^, different from the smaller numbers 
of refs. fraT0|JTl|,|3l. 



^There is a very close link between VMD and the lAM in the vector channels [}liy24| 
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The resulting non-strange and strange normalized scalar form factors of the pion and the kaon are shown in FigJJ, 
Fig.^ and Fig.^. In the case of the non-strange scalar form factor, we show for comparison in Fig^ the one- and 
two-loop CHPT ||l^,Q as well as the dispersion theoretical results and the exponentiated two-loop CHPT 
result Pq] , In the latter case, the two-loop CHPT result is improved by making use of an Omnes resummation in 
terms of the next-to-leading order CHPT phase shifts. Our result is close to the ones obtained by a different method 
in Rcf. and even closer to the exponentiated two-loop CHPT results of rcf. js^. The agreement is worse when 



comparing our results with the so called modified-Omnes representation of refs. }38 13 1. We also remark that the 
two-loop representation covers the main feature of this quantity below W ~ 600 MeV, as it is known since long . 
The strange scalar form factor of the pion is reasonably well described for energies below 350 MeV. In contrast, the 
strange and the non-strange scalar form factor of the kaon are poorly described at one loop, as expected from the 
larger mass of the kaon. 

In Figs.0 and ||, we show the curves from ours fit to the 7r"'"7r~ event distribution in comparison with the experimental 
data from DM2 and MARK-111 for the 25 and the 10 MeV bins, respectively. The data of MARK-HI have been 
multiplied by the factor dm2/^0 mark-iii — l-H'^ in order to facilitate the comparison between both sets of data. 
The agreement with the experimental data is very good as indicated by the low x^/dof of 0.92. By comparing Fig.0 
with the left panel of Fig.|| we see that the event distribution of the two pions is dominated by the strange scalar 
form factor of the pion. Moreover, in FigJ^ the changes in the results when allowing a change in the cut-off qmax by 
±0.1 GeV are shown to be quite small. In Fig|| our prediction for the K^K~ event distribution is depicted. 
Incidentally, we find better agreement with the data of MARK-III than with the ones of DM2 with respect to this 
decay mode. This was also noted in Refs. |l||^. In fact, in these references a fit of similar quality to the data of 
DM2 and MARK-III is also given. The important difference between their method and ours is that we have devised 
a dynamical approach which means that the parameters that enters in our description of the problem are not specific 
to it and can be related to many other physical observables. This is particularly true for and Lg. But even for 
we will see in the next subsection how it can be related to the whole set of U(3) processes that follows from the 
decays of the J/\l/ resonance to any vector resonance belonging to the lightest nonet of vector resonances and two 
pseudoscalars. On the other hand, in ref. HM] no attempt was done to describe the J/^ ^ uitttt decays. 



B. The (jj-meson case 



A priori one can expect that the J/^* wtttt decay requires a rather different dynamical description than that 
for the mode J/^t (pnn considered so far. For instance, the approach of considering the a; as a spectator is by no 
means so clear as for the </> case. Note that the Dalitz plot for this decay has very clear bands due to the decays 
6i(1235) — > UJTT and /2(1270) — > tt+tt" The latter induces a sizeable D-wave contribution so that our approach 

can only be applied for the first few hundred MeV of the two pion event distribution. However, making use of SU(3) 
symmetry, we can extend our considerations from section || and we will present our calculated S-wave contribution 
to the invariant mass distribution of the pions in the J/^t wtttt decay. This calculation is completely fixed in terms 
of the parameters already given for the </) case, cf. Eq.(|44|), except for the global normalization constant for which we 
have also used the experimental data from the J/^ ^ wtt+tt" decays to further constraint its value. 

Let us now be more specific and work out the aforementioned relation between the two cases. An invariant SU(3) 



f 8') -L 

Lagrangian involving the octet and singlet of vector resonances and V^''' , respectively, and the corresponding 
ones of the scalar sources S^*-* and S*^^-* can be written as (making also use of Lorentz invariance) 



.(1) 



(45) 



where g is an overall coupling constant whose precise value is not needed in the following. We only need to determine 
the relative strength of the octet to singlet couplings given in terms of the real parameter v. The symbol (...) refers 
to the trace over the SU(3) indices of the matrices V^^'' and S^^\ These are defined via 







P 

K*- 



P 

1 I It/ 



K* + 



(46) 



and similarly for the S'^^^ matrix. In the previous equation we have denoted by Vs the / = state of the octet of 
vector resonances. This formalism is in close analogy with the one used in CHPT for the octet of pseudoscalars. 
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compare Eq.(^_lJ). Denoting by Sg the 7 = operator of the octet of scalar sources, we can write the terms involving 
Vs and Vi of Eq.dl) as: 

(Vs,aS8 + iyVy^S'^^'>^ . (47) 
Considering ideal mixing between the Vg and the V'^^^ then: 

In an analogous way we will also introduce the scalar sources and defined by 

Note that in a quark model language, consistently with the transformation properties under SU(3), we can write: 

Sij, — ss and S^j — —={uu + dd) . (50) 
v2 

Rewriting Eq.(^7|) in terms of lo, (p, S^^ and S^, we have: 

+ V, + . (51) 

In this way, the parameter introduced in Eqs.(^,^ and fitted in the previous subsection, can now be expressed as: 

„ ^/2(^-l) 

^From this equation we can isolate v — ^'(A^) and then predict the corresponding parameter Aj^, 

A. = -1±M^ (53) 
^/2(^.(A^) - 1) 

We can also obtain from Eq. ( ^T]) the global normalization constant for the oj in terms of the one of the (p since 



Ca 2 + 1/ 



- . (54) 



In Fig.|T^ we show the calculated S-wave contribution to the event distribution of pions in the ^ wtt+tt" decay. 
This calculation does not introduce any new free parameter since all of them have been fixed in terms of the one of 
the (j) case by Eqs.(p3|) and ([54|). The description of the data of MARK-III up to aroimd 0.7 GeV is very good. For 
higher energies the D-wave contribution cannot be further neglected. 



In Ref. the departure from the ideal mixing in the u) — (j) system is thoroughly studied comparing different models, and 
is described by a rotation of the ideal mixing states with a rotation angle \Sv\ ~ 3°. This departure would produce corrections 
of the order of a 5% with respect the ideal mixing situation considered here. However, this value should be compared with the 
departure from 1 of the parameter u which from Eq. ( p2[ ) , taking into account the value of given in Eq. @, is about 40%, 
a much bigger effect than any expected deviation from the ideal mixing situation in the uo — (j) system. Note that v = 1 is the 
value expected from U(3) symmetry and for this case = and A^^ = oo, see Eqs. ( |52|5^ ) 
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V. CONCLUSIONS 



In this work we have addressed the problem of the J/^P decays into a vector lu) and two pseudoscalar mesons 
(Goldstone bosons) measured at DM2 Q and MARK-III These processes are considered to be mediated by 

the corresponding scalar form factors of the pseudoscalar mesons if one considers the emitted vector meson as a 
spectator. Consequently, these reactions are rather interesting since they are very sensitive to OZI violation physics, 
in our scheme parameterized by the constants A^, see Eq.(y), and the low energy constants L4 and Lg of chiral 
perturbation theory. The first of these constants parameterizes the direct admixture of non-strange quarks to the 
scalar interpolating field for our model of the J/\I' decay with the (j) playing the role of a spectator, see FigsJ^,||. The 
two low energy constants enter the one loop description of the pion and kaon scalar form factors. To describe these 
properly for the range of energies relevant here, we have combined information coming from next-to-leading order 
(one loop) chiral perturbation theory (CHPT) with the unitarity requirements which are valid to all orders in the 
chiral expansion. In addition, we also have calculated for the first time the next-to-leading order CHPT kaon scalar 
form factors, for strange, ss, and non-strange, uu -\- dd, scalar-isoscalar quark densities. The unitarity requirements 
were imposed by using the strong / = tttt and KK amplitudes derived in Ref. [Q. The amplitudes given in that 
paper not only describe accurately the S-wave / = and I — 1 strong scattering data but also have been used to 
successfully reproduce or even predict experimental data for the whole set of reactions listed in the Introduction. With 
this input, we have successfully described, from threshold up to around 1.2 GeV,^ the event distribution of the tt+tt" 
system in the J/^* 07r+7r~ decay. We have then predicted, in agreement with the data from MARK-III, the event 
distribution of kaons in the J/^ — + (pK^K^ reaction and the low energy part, where the S-wave dominates, of the 
event distribution of 7r+7r~ pairs in the J/^' — s- wtt'^'tt" decay. Furthermore, the OZI violation parameter comes out 
different from zero. This also holds for the low energy constants L4 and Lg. While the value of the latter agrees with 
previous estimates |l^j3^ , our result for is sizeably larger in magnitude as most previous estimations [ 
However, it is compatible within errors with the quite constraint value derived by combining the information from 
0(/) SU(2) ||l3| and SU(3) ^ CHPT. This offers another indication that the OZI rule does not account for the 
physics in the scalar 0"*'"'' channel, as stressed e.g. in Refs. [^j-^. The scheme employed here offers a unique approach 
to describe the scalar sector, which has been at the heart of many investigations over the last decade. 

Clearly, to further improve the approach presented here, it would be mandatory to not only have event distributions 
but rather normalized data. This would allow one to pin down the low energy constants L4 and Lg more precisely 
together with the OZI violation parameter A^ as well as the product {gruq Bq)^ (i.e. normalization of the scalar form 
factors times the strength of the scalar source to vector meson coupling). With respect to the former, one could also 
try to do a pure one (or even two) loop CHPT calculation for small two-pion invariant masses, i.e. on the left wing 
of the /o(980) resonance. Clearly, the presently available data are not precise enough for successfully doing that, but 
such a program is in principle of the similar interest as the study of chiral dynamics in r decays, see e.g. Refs. p3| , p[ , 
specially when referring to the scalar sector. 
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®In principle, one could also go to higher energies but that is much more difficult due a) to the appearance of multiparticle 
states and b) due to novel interaction vertices as discussed in section F^. 
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FIGURES 




a) b) 



FIG. 1. Quark line diagrams for decay ,7/5' into the (j) and a meson pair (tttt or KK). The quark flavors are explicitly given, 
n refers to the light non-strange u, d quarks. The hatched blob in a) depicts the final state interactions in the coupled m/ KK 
system. 




M 



FIG. 2. Anatomy of the J/^ decay into a 4> and a Goldstone boson pair. S is the interpolating scalar field described in the 
text. The cross-shaded blob symbolizes the final state interactions in the coupled 'kit /KK system. 




a) b) c) d) 



FIG. 3. Feyrunan diagrams for the calculation of the scalar form factors at leading and ncxt-to-loading order in CHPT. The 
scalar source is indicated by the dashed line and the solid lines refer to the psuedoscalars (pions and kaons). From left to 
right: a) Lowest order, b) tadpole contributions, c) unitarity contributions and d) local contact terms from the 0{p*) CHPT 
Lagrangian. The full circles indicates that the vortices come from the lowest order CHPT Lagrangian and the full squares 
symbolize an insertion from next-to-leading order. The wave function reriormalization diagrams are not drawn. 
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FIG. 4. Normalized non-strange scalar form factor of the pion (wide solid line). Dotted and dot-dashed lines: One and two 
loop CHPT results, respectively. The three dashed lines are the dispersion theoretical results from Ref. |39j. The thin solid 
lines represent the exponentiated two-loop CHPT results Our results correspond to the thick solid lines. Left/right panel: 
Real/imaginary part. 
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FIG. 5. Normalized non-strange scalar form factor of the kaon (solid line) in comparison to the one loop CHPT result 
(dashed line). 



18 




FIG. 6. Normalized strange scalar form factor of the pion (left panel) and of the kaon (right panel). Solid lines: Chiral 
unitary approach. Dashed lines: one loop CHPT result. 
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FIG. 7. TT+TT event distribution in the J /'^ —> (pn^-K decay. The width of the bin is 25 MeV. The solid line corresponds 
to the fit Eq.(Q with gmax =0.9 GeV. The dashed line is the best fit with gmax = 1 GeV and analogously the dashed-dotted 
line for gmax =0.8 GeV. 
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FIG. 8. -K^Tv event distribution in the J/^ (p-K^n decay around the /o(980) mass. The width of the bin is 10 MeV. 
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FIG. 9. K event distribution in the J/^ — + (f>K'^ K decay. The upper panel corresponds to the data from DM2 
20 MeV bins. The lower one corresponds to the data from MARK-III Q, 25 MeV bins. 
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FIG. 10. TT~^n event distribution in the J/vt — > ujn'^-K decay. The width of the bin is 25 MeV. Only the S-wave contribution 
is calculated. The onset of the D-wave contribution can be seen at energies larger than 700 MeV. 
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